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Abstract. In many models of dusty objects in space the grains are assumed to be composite or fluffy. However, the computation 
of the optical properties of such particles is still a very difficult problem. We analyze how the increase of grain porosity 
influences basic features of cosmic dust — interstellar extinction, dust temperature, infrared bands and millimeter opacity. 
It is found that an increase of porosity leads to an increase of extinction cross sections at some wavelengths and a decrease at 
others depending on the grain model. However, this behaviour is sufficient to reproduce the extinction curve in the direction of 
the star a Sco using current solar abundances. In the case of the star ( Oph our model requires larger amounts of carbon and 
iron in the dust-phase than is available. Porous grains can reproduce the flat extinction across the 3-8 jim wavelength range 
measured for several lines of sight by ISO and Spitzer. 

Porous grains are generally cooler than compact grains. At the same time, the temperature of very porous grains becomes 
slightly larger in the case of the EMT-Mie calculations in comparison with the results found from the layered-sphere model. The 
layered-sphere model predicts a broadening of infrared bands and a shift of the peak position to larger wavelengths as porosity 
grows. In the case of the EMT-Mie model variations of the feature profile are less significant. It is also shown that the millimeter 
mass absorption coefficients grow as porosity increases with a faster growth occurring for particles with Rayleigh/non-Rayleigh 
inclusions. As a result, for very porous particles the coefficients given by two models can differ by a factor of about 3. 

Key words. Scattering - dust, extinction - comets - interplanetary medium - Stars: individual: £ Oph, a Sco - Galaxy: center 
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c*2 . 1 . Introduction Using the DDA, Voshchinnikov et al. (2005 1 examined the 



ability of the EMT-Mie approach to treat porous particles of 
The steady decrease of estimates of metal abundances in the . . . ™ , ... , .. , 

3 different structure. They show that the latter approach can give 

solar atmosphere over the last years (Holweger 2001 Lodders , .. , . u , .. , ■ 

F 3 v 6 1 1 relatively accurate results only if the very porous particles have 

2003 Asplund et al. 2004i is a serious challenge not only to ,■ ■ ■., , f • ■, . ,. 

' 1 F 1 ' & 3 small (in comparison with the wavelength of incident radia- 

solar physics, but also to dust modelling. This calls for new .. , 1<T , , . , „ . , . , , 

v 3 b tion) Rayleigh inclusions. Otherwise, the approach becomes 

dust models able to produce the same extinction with a smaller . , .. , n c . . 

r inaccurate when the porosity exceeds -0.5. At the same time, 
amount of solid material. A solution to the problem could be .. , .. . , . , .. , , 

v the optical properties of heterogeneous spherical particles hav- 

provided by an "admixture of vacuum" i.e. by the porosity of • • , • t ■ ■ /TJ , • , ■ D , • , ^ 

f 3 ing inclusions of various sizes (Rayleigh and non-Rayleigh) 

in ers e ar grains. anc j large porosity were found to closely resemble those 

Grain aggregates with large voids can form during the of spheref , wkh & largg number ( > 15 _ 2Q) Qf different ky _ 

growth of interstellar grains due to their coagulation in dense erg The errQrs - m extinction efficiency factors are smaller than 

molecular cloud cores (Dorschner & Henning [1995}. The in- 1Q _ 2Q% tf ±e gize parameter < 15 and porosity is equal to 

ternal structure of such composite grains can be very compli- Q g Note ^ Ms consideration was restricted by spheres, not 

cated, but their optical properties are often calculated using the yery absorbing materials (silicate and amorphous cai'bon) and 

Mie theory for homogeneous spheres with an average refractive ^ scattering characteristics (extinction, scattering, ab- 

index derived frorneflfective medium theory (EMTjsee, e.g., ion efficie factor albedo and asymme t ry parameter) 

Mathis & Whiffen[I»ones[l988l OssenkopfEHl Mathis but P not the differe y ntial crQSS sectionfj or e y lements y o F f the S J_ 

1996 Li & Greenberg 1998 Win & Krivova|2000). tering matrix NeverthelesSj very simple comput ational models 

Another approach to calculate the optical properties of such instead of time -consuming DDA calculations give us a useful 

aggregates is the application of complex, computationally time way to ^ composite gra ins of different structure, 
consuming methods such as the discrete dipole approximation 

(DDA; see, e.g., WrightQ987 Kozasa et al. 1992"! 1 19931 Wolff In this paper, we apply the particle models of porous inter- 

et al. 1 19941 Stognienko et al. 1 19951 Kimura & Mann 1998 1. stellar dust grains based on the EMT-Mie and layered-sphere 
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calculations. The models described in Sect. |2 are assumed 
to represent composite particles with small (Rayleigh) inclu- 
sions and inclusions of different sizes (Rayleigh and non- 
Rayleigh). The wavelength dependence of extinction is dis- 
cussed in Sect. 13.1 1 Sections |3 .21 and l3~3l contain the applica- 
tion of the models to calculations of the extinction curves in 
the directions of two stars, using new solar abundances and 
the near-infrared (IR) extinction in the directions along the 
Galactic plane. The next sections deal with grain temperatures 
(Sect. 14. 1 1 , profiles of IR silicate bands (Sect. 14.21 . and grain 
opacities at A = 1 mm (Sect. l4~3l . These quantities are espe- 
cially important for the analysis of observations of protoplane- 
tary discs (Henning et al.[2p05 l. Concluding remarks are pre- 
sented in Sect. [5] 

2. Particle models 

Information about the structure of grains can be included in 
light scattering calculations directly (layered spheres) or can 
be used to find the optical constants (EMT-Mie model). We 
consider models of both types. 

Following previous papers (Voshchinnikov & Mathis 1999 
Voshchinnikov et al. 2005 1, we construct layered grains as 
particles consisting of many concentric spherical layers of 
various materials, each with a specified volume fraction V, 
(£,-V,y V tota i = 1). Vacuum can be one of the materials, so a 
composite particle may have a central cavity or voids in the 
form of concentric layers. The number of layers is taken to be 
1 8 since Voshchinnikov et al. ( 2005 1 have shown that this was 
enough to preclude an influence of the order of materials on 
the results. For a larger number of layers, one can speak of the 
optical characteristics determined by the volume fractions of 
different constituents only. 

In the case of the EMT-Mie model, an average (effec- 
tive) refractive index is calculated using the popular rule of 
Bruggeman (see, e.g., Chylek et al. 2000] Kriigel I2003> . In 
this case, the average dielectric permittivity e e ff (1 ' is calculated 
from 

i 

where = V,7Vt tai is the volume fraction of the ith material 
with the permittivity e, . 

The amount of vacuum in a particle can be characterized by 
its porosity f (0 < f < 1) 

P = Vvac/Vtota, = 1 - VsoHd/W (2) 

To compare the optical properties of porous and compact parti- 
cles, one should consider the porous particles of radius (or size 
parameter) 

'"compact '"compact 
^ous = (1 _ ry/3 = (ysol . d/ y toml) l/3- W 

As "basic" constituents, we choose amorphous carbon 
(AC1; Rouleau & Martin I199U and astronomical silicate 

1 The dielectric permittivity is related to the refractive index as s = 

2 

111'. 



(astrosil; Laor & Draine 1993 1. The refractive indices for 
these materials and some others considered in Sect. [5] were 
taken from the Jena-Petersburg Database of Optical Constants 
(JPDOC) described by Henning et al. (1999) and Jager et 
al. 

The application of the standard EMT is known to corre- 
spond to the case of particles having small randomly located in- 
clusions of different materials (Bohren & Huffman 19831. The 
optical properties of such particles have been well studied (see, 
e.g., Voshchinnikov 2004 and references therein). 

However, one needs the DDA or other computationally 
time consuming techniques to treat particles with inclusions 
larger than the wavelength of incident radiation. The difference 
in the optical characteristics of particles with small and large 
inclusions has been discussed in previous studies (e.g., Wolff et 
al. fT994l fT998l. The fact that this difference drastically grows 
with the porosity f and becomes quite essential already for 
V ^ 0.5 has been discovered only recently by Voshchinnikov 
et al. (2005). They also found that the scattering properties of 
particles with inclusions of different sizes (including those with 
sizes comparable to the wavelength) were very close to those 
of layered particles, having the same size and material volume 
fractions. A similar conclusion was reached for an ensemble of 
particles where each particle has inclusions of one (but differ- 
ent) size only. In both cases of "internal" or "external" mixing, 
the model of layered spheres can be applied. 

These results are of particular importance for applications 
including the astronomical ones where we deal with very 
porous particles. Note that instead of time consuming calcu- 
lations with DDA-like codes, one can use the exact solution 
to the light scattering problem for layered spheres, which is as 
quick as Mie theory, to estimate the properties of particles with 
Rayleigh/non-Rayleigh inclusions. 

Thus, the models of the homogeneous sphere (with EMT) 
and layered spheres, both having fast implementations, allow 
one to probe the difference in optical properties caused by the 
different structure of scatterers. 



3. Interstellar extinction and interstellar 
abundances 

3.1. Wavelength dependence of extinction 

As it is well known, the wavelength dependence of interstellar 
extinction A{A) is determined by the wavelength dependence 
of the extinction efficiencies Q ex t{A). This quantity is shown 
in Fig. ^foi" particles of the same mass, but different poros- 
ity. The volume fractions of AC1 and astrosil are equal, i.e. 

V AC l/VtotaJ = Vastrosil/Vtotal = 1 /2 (V so lid/ Vtotal) = 1/2(1 ~P). 

The radius of compact grains is r s compact = 0.1 //m. The de- 
pendence of <2 ext on A for compact particles is close to that 
of the average interstellar extinction curve in the visible-near 
UV (1 fim~ l < A < 3/iirT 1 ) where it can be approxi- 
mated by the power law A{A) cc A' 1 33 (see discussion in 
Voshchinnikov [2004}. From Fig. [2 one can conclude that the 
extinction produced by particles with small inclusions and lay- 
ers differs little for compact and slightly porous particles. The 
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difference becomes most pronounced in the near and far-UV 
(for,*- 1 ^2.5-3/mT 1 ). 




Fig. 1. Wavelength dependence of the extinction efficiency fac- 
tor for spherical particles with r s< compact = 0.1 //m. The particles 
are of the same mass but of different porosity. Upper panel: cal- 
culations based on the EMT-Mie theory. Lower panel: calcula- 
tions based on the layered-sphere theory. 
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A , /im 1 

Fig. 2. The same as in Fig.^but now for normalized extinction 
cross section. 



gain in dust-phase abundances and would contradict the wave- 
length behaviour. 

The role of porosity in extinction is better seen from Fig. 13 
where we give the wavelength dependence of the normalized 
cross section 

c („) Cext(porous grain) 

ext C ex t(compact grain of same mass) 



As follows from Fig.^ the wavelength dependence of ex- 
tinction flattens as porosity increases. It is well known (see, 
e.g., Greenberg 1978 1 that different particles produce compara- 
ble extinction if the products of their size r and refractive index 
are close, i.e. 



-2/3 



2ext(porous grain) 



(5) 



r\m — II 



const. 



The average refractive index of particles with a larger fraction 
of vacuum is closer to 1. Despite a larger radius (e.g., from 
Eq. follows that r s = 0.22 /mi if P — 0.9 and r S; compact = 
0.1 yum), the product given by Eq. (0} decreases because of an 
drop of \m — 1|. This implies that a steeper extinction with the 
wavelength dependence closer to oc A -133 is produced by for 
compact particles with larger radii and, consequently, a larger 
amount of solid material. This behaviour is also observed for 
particles of other masses (i.e., compact spheres of other radii). 
Therefore, an interpretation of the observed interstellar extinc- 
tion curve using only very porous grains should not give any 



2ext(compact grain of same mass) 

This quantity shows how porosity influences the extinction 
cross section. As follows from Fig. [2] both models predict a 
growth of extinction of porous particles in the far-UV and a de- 
crease in the visual-near-UV part with growing P. However, 
(4) the wavelength interval where C^J < 1 is narrower and the 
minimum is less deep in the case of layered spheres. In com- 
parison with compact grains and particles with Rayleigh in- 
clusions, particles with Rayleigh/non-Rayleigh inclusions can 
also produce rather large extinction in the near-IR part of spec- 
trum. This is especially important for the explanation of the 
flat extinction across the 3-8 //m wavelength range measured 
for several lines of sight (see Sect. 13 .31 . At the same time, as 
follows from Fig. [3] for particles with r SiCompact > 0.1 /mi the 
normalized UV cross sections grow faster for the Bruggeman- 
Mie theory than for layered spheres. Thus, an addition of vac- 
uum into particles does not mean an automatic growth of ex- 
tinction at all wavelengths and a saving in terms of solid-phase 
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,(n) 
-"ext 



(n) 



r aiCompact =0.0 1/j.m 

r a , [ , ompa( , t =0.0^/;.m 

I'a, compact — 0.05/i.m 

rg,compact = 0-10jU.m 

***** r SiCompaet =0.^0^tm x 

- — — r HiCompact — 0.30,Ltm ^ 



Brugge man— Mie, 5P=0.9 

— ~ \ 

* * + * 



9 10 
X , ^trn 1 



layered sphere, ^=0.9 



X , 

Fig. 3. Wavelength dependence of the normalized extinction 
cross section for spherical particles with the same porosity 
P = 0.9. Upper panel: calculations based on the EMT-Mie 
theory. Lower panel: calculations based on the layered-sphere 
theory. 



elements. Evidently, the final decision can be made after fitting 
the theoretical calculations with observations at many wave- 
lengths. 



determined by Sofia & Meyer (2001 ) for F and G stars must 
be revised downward as it has been done recently for the Sun. 
This should lead to the agreement between abundances found 
for stars of different types. Note also that the current solar abun- 
dances of oxygen and iron (see Table[5} are close to those found 
from high-resolution X-ray spectroscopy. Juett et al. ( 2004 1 in- 
vestigated the oxygen K-shell interstellar absorption edge in 
seven X-ray binaries and evaluated the total O abundances. 
These abundances lie between 467 ppm (2) and 492 ppm. Schulz 
et al. (2002) evaluated the total abundance of iron towards the 
object Cyg X-l to be [Fe/H] cosm i c « 25 ppm. 

We applied the model of multi-layered porous particles to 
explain the absolute extinction in the direction to the two stars. 
A first estimate has been made in order to find a possibility 
to enlarge the extinction per unit mass and to minimize the 
amount of solid phase material. Several materials as compo- 
nents of composite grains were considered. Among the car- 
bonaceous species, the amorphous carbon in the form of Bel 
(Rouleau & Martin n"991t was found to produce the largest ex- 
tinction. Also the extinction of iron oxides strongly increases 
with the growth of porosity. Although there are no very good 
constraints on the abundance of oxides, FeO is considered as a 
possible carrier of the 21 fj.m emission observed in the spectra 
of post-AGB stars (Posch et al. 2004). Very likely, such parti- 



cles particles can be produced in redox reactions (Duley 1980 
Jones 1993. 



A(A) 



■ * ■ - observations 
o o o □ theory 

— — Bel/pyroxene (I) 

graphite (II) 

- - magnetite (III) 
forsterite (IV) 

■ ~ CCM fitting (R,=3.09) 



t Oph 




3.2. Extinction in the directions to the stars £ Oph and 

cr SCO 

The basic requirement for any model of interstellar dust is the 
explanation of the observed extinction law along with the dust- 
phase abundances of elements in the interstellar medium. These 
abundances are obtained as the difference between the cosmic 
reference abundances and the observed gas-phase abundances. 
However, the cosmic abundances are not yet conclusively es- 
tablished and usually this causes a problem. For many years, 
the solar abundances were used as the reference ones, until the 
photospheres of the early-type stars were found not to be as 
rich in heavy elements as the solar photosphere was (Snow & 
Witt 1996 ). These stellar abundances caused the so-called "car- 
bon crisis". Abundances of the most important dust-forming 
elements (C, O, Mg, Si, Fe) required by the dust models were 
larger than available. However, during the past several years 
the solar abundances dropped and now they approach the stel- 
lar ones (see Asplund et al. 2004). Evidently, some abundances 



Fig. 4. Observed and calculated extinction in the direction to £ 
Oph. The errors of the observations are the result of a parame- 
terization of the observations (see Fitzpatrick & Massa 1990). 
The contribution to the theoretical extinction from different 
components is also shown. The dot-dashed curve is the approx- 
imation with the observed value Ry as suggested by Cardelli et 
al. fT989l 



We fitted the observed extinction toward the stars £ Oph and 
cr Sco. For these stars there exist well determined extinction 
curves and gas-phase abundances. It is also important that the 
major part of extinction in these directions is produced in one 
diffuse interstellar cloud (Savage & Sembach 1996, Zubko et 
al. 1996 ). This allows us to exclude possible large variations in 
dust composition along the line of sight. 

2 parts per million hydrogen atoms 
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X , fim 1 

Fig. 5. The same as in Fig.@]but now for cr Sco. The observa- 
tional data were taken from Wegner ( 2002 1. 



Observed and calculated extinction curves are plotted in 
Fig s -E](<r Oph) and|5](cr Sco). As follows from the previous 
Section, the use of only porous or only compact grains appar- 
ently does not result significant benefit in the solid-phase abun- 
dances. Therefore, our models are the combination of compact 
and porous particles. They consist of three or four grain popu- 
lations: 

(I) . Porous composite (multi-layered) particles (Bel — 
5%, pyroxene, Feo.sMgo.sSiC^ — 5% for ( Oph or forsterite, 
Mg 2 Si0 4 — 5% for cr Sco and vacuum — 90%) with 
the power-law size distribution having an exponential decay 
Wd( r s) <* r~ 2 - 5 exp(-10/r s ). The lower/upper cut-off in the 
size distribution is 0.015 fim/0.25 fim and 0.05 //m/0.50 fim for 
( Oph and cr Sco, respectively. 

(II) . Small compact graphite 3 grains with a narrow power- 
law size distribution (ndC^s) <* r 7 2 ' 5 , r s = 0.01 - 0.02 fim). 

(III) . Porous composite grains of magnetite (Fe304 — 2%, 
vacuum — 98% for f Oph and Fe304 — 8%, vacuum — 92% 
for cr Sco) with a power-law size distribution «d( r s) K r 7 2 ' 5 - The 
lower/upper cut-off in the size distribution is 0.005 //m/0.25 fim 
and 0.05 fim/0.35 fim for £ Oph and cr Sco, respectively. 

(IV) . Compact grains of forsterite (Mg2Si04) with the 
power-law size distribution (only for £ Oph, «d('"s) x 7\T 3 5 , 
r s ,min = 0.10 fim, r Sjtnax = 0.25 fim). 

Figures |4] and [5] also contain the extinction curves cal- 
culated using the approximation suggested by Cardelli et al. 
(1989l with the coefficients revised by O'Donnell J1994i . 
Cardelli et al. ( 1989 1 found that the extinction curves from the 
UV through the IR could be characterized as a one-parameter 
family dependent on the ratio of the total extinction to the se- 
lective one Ry = Ay/E(B - V). We used the observed values 
of Ry in order to plot the CCM approximation. It is seen that 
this relation describes quite well the extinction for f Oph but 
not for cr Sco. 



3 The calculations for graphite were made in the "2/3-1/3" approx- 
imation: <2 cxl = 2/3 <2cxt( £ ±) +1/3 <2cxt( £ nX where e ± and e\\ are the 
dielectric functions for two cases of orientation of the electric field 
relative to the basal plane of graphite. 



The contributions from different components to the cal- 
culated extinction are given in Tables ^ an d HI an d shown in 
Figs. 0] and [5] The Tables contain also the dust-phase abun- 
dances of five dust-forming elements for several grain popula- 
tions. They were calculated for ratios of the extinction cross- 
section to particle volume averaged over grain size distribution 
(see Eq. (3.36) in Voshchinnikov 2004). 

Table [3] gives the current solar abundances of five dust- 
forming elements according to Asplund et al. (2004) as well as 
the "observed" (solar minus gaseous) and model abundances. 

The dust-phase abundances in the line of sight to the star 
f Oph (HD 149757) were taken from Table 2 of Snow & 
Witt J1996i . In our calculations, we adopted the following 
quantities for ( Oph: a total extinction Ay = 0^94 (4) , colour 
excess £(B - V) = (J?32 and total hydrogen column density 
N(H) = 1.35 10 2I cnr 2 (Savage & Sembach. [T996l . The ex- 
tinction curve was reproduced according to the parameteriza- 
tion of Fitzpatrick & Massa ( 1990). 

For cr Sco (HD 147165), we used the extinction curve, 
the colour excess E(B - V) = 0^35 and the total extinction 
Ay = 1™13 according to Wegner (2002). The hydrogen col- 
umn density N(H) = 2.46 10 21 crrT 2 was adopted from Zubko 
et al. ( 1996 ). The gas-phase abundances were taken from Allen 
et al. (fT990l . 

The dust-phase abundances required by the model are 
larger than the observed ones in the direction to f Oph (for 
C and Fe) and smaller than the observed abundances in the di- 
rection to cr Sco. Note that for cr Sco the required amount of 
C and Si in dust grains is the lowest in comparison with previ- 
ous modelling. This is due to the use of highly porous particles 
which give considerable extinction in the UV and near-IR (see 
Figs.[2]and[3} and allow one to "save" material. For example, 
the extinction model of cr Sco with compact grains presented 
by Zubko et al. ( fT996l requires 240 - 260 ppm of C and 20 - 
30 ppm of Si and the model of Clayton et al. ( 2003 ) needs 155 
ppm of C and 36 ppm of Si (cf. 137 ppm and 8.8 ppm from 
Tabled}. 

The models presented above are based on the light scatter- 
ing calculations for particles with Rayleigh/non-Rayleigh in- 
clusions (layered spheres). It is evident that the observed ex- 
tinction can be also reproduced if we use the particles with 
Rayleigh inclusions (i.e., if we apply the EMT-Mie theory). 
Our estimates show that despite a larger extinction in the UV 
this model requires more material in the solid phase in compar- 
ison with layered spheres because of a smaller extinction in the 
visual-near IR part of the spectrum. 

3.3. Near infrared extinction in the Galactic plane 

We now consider the possibility of explaining the flat extinc- 
tion across the 3-8 fim wavelength range observed for sev- 
eral lines of sight. This flattening was first measured by Lutz 
et al. (1996) toward the Galactic center with ISO, using hy- 
drogen recombination lines. Later Lutz Q999) confirmed the 

4 This value was obtained from the relation A v = 1.12 E(V - K) 
(Voshchinnikov & irin[T987l and a colour excess £(V - K) = (T84 
(Serkowski et al.[T575J. 
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Table 1. Contribution of different grain populations to Ay and dust-phase abundances for the model of f Oph (in ppm) 



Component 




c 


O 


Mg 


Si 


Fe 


(I) Bel/pyroxene/vacuum 


a; 1 5 8 


123 


68 


11.3 


22.6 


11.3 


(II) Graphite 


0P075 


96 










(III) Magnetite/vacuum 


a? 22 




33 






24.8 


(IV) Forsterite 


(F065 




23 


11.4 


5.7 




Total 


a? 94 


219 


124 


22.7 


28.2 


36.1 



Table 2. Contribution of different grain populations to Ay and dust-phase abundances for the model of <x Sco (in ppm) 



Component 


A v 


C 


O 


Mg 


Si Fe 


(I) Bel/forsterite/vacuum 


(T50 


58 


35.4 


17.7 


8.8 


(II) Graphite 


Of 11 


79 








(III) Magnetite/vacuum 


(T52 




35.4 




26.6 


Total 


r?13 


137 


71 


17.7 


8.8 26.6 



Table 3. Observed and model dust-phase abundances (in ppm) 



Element 



Solar* 
abundance 



obs model 



cr Sco 
obs model 



c 


245 


110 


219 


176 


137 


o 


457 


126 


124 


85 


71 


Mg" 


33.9 


31.9 


22.7 


30.9 


17.7 


Si 


34.2 


32.6 


28.2 


32.4 


8.8 


Fe 


28.2 


28.2 


36.1 


27.9 


26.6 



* According to Asplund et al. 1 2004 1. 

** The abundance of Mg was recalculated with the oscillator strengths 
from Fitzpatrick (I1997> . 



effect using more recombination lines. Recently, Indebetouw 
et al. (2005 1 found a similar flat extinction along two lines of 
sight: / = 42° and I - 284°. The extinction was obtained at 
seven wavelengths (1 .2 - 8 /mi) by combining images from the 
Spitzer Space Telescope with the 2MASS point-source catalog. 

We used the average extinction given in Table 1 of 
Indebetouw et al. (2005 1 and transformed it into magnitudes 
of extinction per kpc using the measured value — 07T5 + 

0?1 kpc -1 (Indebetouw et al. |2005| . The observations are plot- 
ted in Fig. [6] together with three theoretical curves. Because 
we have little information about the UV-visual extinction and 
gas-phase abundances in these directions, we (rather arbitrar- 
ily) applied the model used for £ Oph (porous component (I): 
Bel/pyroxene, porosity 90%; see Sect. 13. 2\ . This model (short 
dashed curve on Fig. |6j well explains the flat extinction at 
A > 3/mi (5) but the extinction in the J and H bands is too 
small. Compact particles (long-dashed curve in Fig. [6} produce 
even larger extinction at these bands than the observed one. 



5 Note that porous particles from magnetite (component (III); see 
Sect. 13.21 cannot fit well the extinction at these wavelengths because 
of a bump at A « 2 jim. 



O 



0.6 



0.5 



iWWrfVw- observations 

,^ = 0.0 

3> = Q 

^ = 0.9 

• CCM fitting (R v =3.1) 



0.4 
R 0.3 



P 



0.2 



< 0.1 
0.0 




0.0 0.2 



0.4 



0.6 0.8 



A" 



1.0 



Fig. 6. Observed and calculated extinction in the near-IR part of 
spectrum. The observations correspond to the average extinc- 
tion for two lines of sight along the Galactic plane (Indebetouw 
et al. 2005 1 transformed into magnitudes of extinction per kpc. 
The theoretical extinction was calculated for component (I) of 
the model used for £ Oph CP = 0.9, short dashed curve). Two 
other curves correspond to the same component but with an- 
other particle porosity. The dot-dashed curve is the approxima- 
tion of Cardelli et al. dT9"59l with R Y = 3.1. 



However, the extinction from such particles at longer wave- 
lengths decreases rapidly. Our preliminary analysis shows that 
particles with a porosity of about 0.6 (solid curve on Fig. |6j 
can be chosen as an appropriate model. Evidently, a similar 
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curve can be obtained as a combination of compact and very 
porous particles. Quite close extinction can be found with the 
CCM approximation and the standard value Ry = 3.1. We ar- 
bitrarily extrapolated this approximation to long wavelengths 
(A~ l < 0.3 /UirT 1 ) where it gives too small extinction. 

Extinction produced by porous grains was also rather flat 
between 1.0 and 2.2 fim (for example, A(A) oc A~ 13 for P = 
0.6) as was detected for several ultracompact HII regions with 
Ay £ 15 m (Moore etal.[ 



In order to calculate the dust-phase abundances for mod- 
els presented in Fig. [6] we estimated the total hydrogen 
column density using Eqs. (3.26), (3.27) and (3.22) from 
Voshchinnikov (2004i. First, we found the column density of 
atomic hydrogen A^(HI) from the extinction at J band and then 
transformed A^(HI) into a total hydrogen column density N(H) 
using the ratio of total to selective extinction Ry — 3 . 1 . A value 
of N(H)/D = 2.79 10 21 cm^kpc" 1 was obtained. The calcu- 
lated abundances are given in Table 0] which also contains 
the visual extinction calculated for three models. Note that the 
model of grains with porosity P = 0.6 gives the largest contri- 
bution to A v in comparison with two other models. 



Table 4. Dust-phase abundances (in ppm) and visual extinction 
for three models presented in Fig.|6] 



Element 


P = 


P = 0.6 


IP = 0.9 


C 


96 


84 


62 





53 


46 


34 


Mg 


8.9 


7.7 


5.7 


Si 


18 


16 


11 


Fc 


8.9 


7.7 


5.7 


A v 


OP60 


0P79 


W61 



M 
CD 



CD 



CD 

F— I 

a 

CD 



Bruggeman— Mie 

3>=0.0 

o o o oo !P=0. 1 
3>=0.3 

MHHM 5>=0.5 

***** 5>=0.7 
^=0.9 




0.01 




0.1 

r s, compact! /Xtn 


layered spheres 










/ ° li 
o 






. • «, 














>f * * * 


* 


* * 


0.01 




0.1 



^s, compact* ^Tfl 

Fig. 7. Size dependence of the temperature for spherical parti- 
cles. The particles are located at a distance of 10 4 /?* from a 
star with an effective temperature = 2500 K. Upper panel: 
calculations based on the EMT-Mie theory. Lower panel: cal- 
culations based on the layered-sphere theory. 



There are useful observational data for foreground stars 
in the field / = 284 °. For stars HD 90273, HD 93205 and 
HD 93222, Wegner ( 120021 ) and Barbara et al. (120041 esti- 
mate the values of Ry which lie between 3.4 and 4.0. For 
HD 90273 Barbara et al. (2004) also find an anomalously high 
gas to dust ratio JV(H)/£(B-V) = 1.04 10 22 atoms cm" 2 mag." 1 
Enlargement of Ry increases the required dust-phase abun- 
dances while the decrease of the gas to dust ratio reduces them. 
Andre et al. (2003 ) measured the interstellar gas-phase oxygen 
abundances along the sight lines toward 5 early-type stars with 
/ = 285°3 - 287°7 and b = -5°5 - + 0°1. The values of 
[0/H] g = 443 ppm vary from 356 ppm to 512 ppm, the aver- 
age value being 443 ppm. This gives for the mean dust-phase 
abundance [0/H]d = 14 ppm. However, the extinction curves 
for HD 93205 and HD 93222 published by Wegner J2002l l 
have strong UV bumps and flat extinction in the far-UV. This 
means that extinction can be mainly produced by carbonaceous 
grains. Evidently, the most reasonable way to solve the prob- 
lem of abundances is a re-examination of the reference cosmic 
abnudances and a detailed study of their local values. 



4. Infrared radiation 

4.1. Dust temperature 

The commonly used equilibrium temperature of cosmic grains 
is derived from a balance between the energy gain due to ab- 
sorption of the UV and visual stellar photons and the energy 
loss due to re-emission of IR photons. The temperature of 
porous and compact particles of different size and porosity is 
shown in Figs. and [8] as a function of particle size and poros- 
ity, respectively. The results were calculated for particles lo- 
cated at a distance of 10 4 /?* from a star with an effective tem- 
perature T+ = 2500 K. In the case of the layered spheres, an 
increase of the vacuum fraction causes a decrease of the grain 
temperature if the amount of the solid material is kept constant. 
This behaviour holds for particles of all sizes as well as for par- 
ticles located closer to the star or farther away and for other val- 
ues of T+. If the EMT-Mie theory is applied, the temperature 
drops when the porosity grows up to ~ 0.7 and then starts to 
increase (see Fig.0 upper panel and Fig. [8}- Such a behaviour 
corresponds to the results of Greenberg & Hage (1991) who 
found an increase of temperature for large grain porosity (see 
Fig. 4 in their paper). 
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0.0 0.2 0.4 0.6 0.8 1.0 

porosity 



Fig. 8. Dependence of the dust temperature on particle porosity. 
Solid lines: calculations based on the EMT-Mie theory, dashed 
lines: calculations based on the layered-sphere theory. Other 
parameters are the same as in Fig.0 



As it follows from Figs.0and[8] the difference in the tem- 
perature of very porous grains calculated using the two models 
can reach ~ 6 K or ~ 15% while the temperature of compact 
CP = 0) composite grains differs by less than 1%. Note that the 
relative difference in temperatures of ~ 15% between particles 
of the two types is kept for other stellar temperatures (e.g., in 
the case of the Sun or an interstellar radiation field). 

The intermediate porosity of grains leads to a shift of the 
peak position of IR emission to larger wavelengths in compar- 
ison with compact particles. This occurs independently of par- 
ticle structure (small or various size inclusions). However, very 
porous grains with Rayleigh/non-Rayleigh inclusions are ex- 
pected to be systematically cooler than particles with Rayleigh 
inclusions. Such a difference can be of great importance 
at a lower temperature regime because it can influence the 
growth/destruction of mantles on grains in molecular clouds. 

4.2. Infrared features 

It is well known that the shape of the IR dust features is a good 
indicator of the particle size and chemical composition. With 
an increase of the size, a feature becomes wider and eventu- 
ally fades away. For example, in the case of compact spheri- 
cal grains of astrosil, the 10 /mi and 18 /mi features disappear 
when the grain radius exceeds ~ 2 - 3 /mi. Observed differ- 
ences in small scale structure of the features are usually at- 
tributed to variations of the composition (e.g., changes of the 
ratio of magnesium to iron in silicates) or material state (amor- 
phous/crystalline). 

In Fig. [9] we compare the wavelength dependence of the 
absorption efficiency factors for particles of the same mass but 
different structure. The upper panel shows results obtained with 
the EMT-Mie model for particles with Rayleigh inclusions. 
It can be seen that the central position and the width of the 
dust features does not really change. Larger changes occur for 
the layered-sphere model (Fig. |9j lower panel). In this case a 
growth of P causes a shift of the center of the feature to longer 
wavelengths and its broadening. For particles with P = 0.9, the 
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Fig. 9. Wavelength dependence of the absorption efficiency fac- 
tors for spherical particles of radius r Sj0ompact = 0.1 /mi. Upper 
panel: calculations based on the EMT-Mie theory. Lower panel: 
calculations based on the layered-sphere theory. 



10 /mi feature transforms into a plateau while the 18 /mi feature 
disappears. 

We plotted in Fig. EI|our data from Fig. |5]in a normal- 
ized manner together with observations of (3 Pictoris made by 
Knacke et al. JT9931 and Telesco & Knacke fTWTY As follows 
from Fig.^| for given optical constants of the silicate the ob- 
served shape of the 10 //m feature is better reproduced by either 
compact or porous particles with small size inclusions of mate- 
rials. Note that in the case of /? Pictoris the EMT-Mie calcula- 
tions were earlier used by Li & Greenberg ( 1998 1 for the expla- 
nation of the 10 /mi emission feature and by Voshchinnikov & 
Kriigel ( 1999 1 for the interpretation of the positional and wave- 
length dependence of polarization. The best fit was obtained for 
very porous particles: P « 0.95 and 0.76, respectively. 

Figure shows the normalized absorption efficiency fac- 
tors for spherical particles of radius r SjCompact = 2 /mi. For com- 
pact grains the 10 /mi feature almost disappears. When the 
porosity increases the strength of the feature grows in the case 
of the Bruggeman-Mie calculations. This tendency coincides 
with the results shown in Fig. 7 of Hage & Greenberg ( 1990l 
who found that the higher the porosity, the sharper the silicate 
emission became. For the case of layered spheres the feature 
becomes only slightly stronger but its peak shifts to longer 
wavelengths. 
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Fig. 10. Emission in the disc around the star B Pictoris in the re- 
gion of silicate 10 pm band. Stars and squares are the observa- 
tions of Knacke et al. < 



( 1993 ) and of Telesco & Knacke ( 1991 1. 
The curves present the results of calculations for particles of 
radius r s>comp act = 0.1 pm as shown in Fig.|9]but normalized at 
A = 9.6 pm. 



A standard "compact" approach to the modelling of the 
10 pm feature was used by van Boekel et al. ( 20031 12005> 
and Przygodda et al. (|2003 ) who considered the flattening of 
the 10 pm feature as an evidence of grain growth in the discs 
around Herbig Ae/Be stars and T Tauri stars, respectively. Our 
investigations show that the variations of the shape of the fea- 
ture and its position and strength can also be attributed to the 
change of porosity and relative amount of carbon in composite 
grains of small sizes. 



4.3. Dust opacities 

The dust opacity or the mass absorption coefficient of a grain 
material k(A) enters directly in the expression for the dust mass 
of an object M& which is determined from optically thin mil- 
limeter emission 



M d = 



F mm (A)D 2 
K(A)B A (T d y 



(6) 



Here, F mm (A) is the observed flux, D the distance to the ob- 
ject, B^iTi) the Planck function, 7d the dust temperature. The 
mass absorption coefficient k(A) depends on the particle volume 



Fig. 11. Wavelength dependence of the normalized absorption 
efficiency factors for spherical particles of radius r s , compact = 
2 pm. Upper panel: calculations based on the EMT-Mie theory. 
Lower panel: calculations based on the layered-sphere theory. 



Vtotai, m e material density pa and the extinction cross-section 
Cext as follows: 



K(A) = 



Cp- 



Pd^total 



3_ 
Pi 



A j I e efF + 2 



(7) 



At long wavelengths the scattering can be neglected (C ex t ~ 
C a bs) and C a bs can be evaluated in the Rayleigh approximation. 
Then the mass absorption coefficient does not depend on the 
particle size as shown in the right part of Eq. Q. The effec- 
tive dielectric permittivity e e $ in Eq. Q can be found from the 
Bruggeman rule (see Eq. (Q)) or the layered-sphere rule of the 
EMT (see Eqs. (7), (8) in Voshchinnikov et al. l2005l . 

Extensive studies of the mass absorption coefficient depen- 
dence on the material properties and grain shape are summa- 
rized by Henning 1 1996 1 who, in particular, notes that the opac- 
ities at 1 mm are considerably larger for non-spherical particles 
than for spheres (see also Ossenkopf & Henning ll9941 l. We find 
that a similar effect (an increase of opacity in comparison with 
compact spheres) is produced by inclusion of a large fraction of 
vacuum into the particles. This follows from Table[5]where the 
opacities at A = 1 mm are presented. This Table contains the 
results for particles consisting of three materials (AC1, astrosil 
and vacuum) or two materials (AC1 or astrosil and vacuum). In 
the first case, the volume fractions of AC1 and astrosil are equal 
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Table 5. Mass absorption coefficients at A = 1 mm of compact and porous spheres consisting of AC1* and (or) astrosil" 



p 




AC1 + astrosil 




astrosil 






AC1 




Pi 


k, cm 2 /g 


Pd 


k, cm 2 /g 


Pi 


K, 


cm 2 /g 




Brugg.-Mie lay. spheres 




Brugg.-Mie 


lay. spheres 




Brugg.-Mie lay. spheres 


0.00 


2.58 


1.58 1.58 


3.30 


0.310 


0.310 


1.85 


4.37 


4.37 


0.10 


2.32 


1.89 1.55 


2.97 


0.371 


0.334 


1.66 


5.13 


4.60 


0.30 


1.80 


2.75 1.87 


2.31 


0.548 


0.446 


1.30 


7.09 


5.77 


0.50 


1.29 


3.83 2.57 


1.65 


0.778 


0.646 


0.925 


9.22 


7.88 


0.70 


0.772 


3.94 4.04 


0.990 


0.794 


1.05 


0.555 


9.14 


11.9 


0.90 


0.258 


2.45 8.12 


0.330 


0.431 


2.20 


0.185 


5.94 


21.8 



m(A = 1mm) = 2.93 + 0.276/, p d = 1.85 g/cm 3 
* m(A = 1mm) = 3.43 + 0.050/, p d = 3.3 g/cm 3 



(VAci/V to tai = Vastrosii/Vtotd = 1/2(1- <P)) while in the second 
case the volume fraction of solid material is 1 - V. It can be 
seen that the values of k are generally larger for particles with a 
larger fraction of vacuum. This is related to the decrease of the 
particle density pd which is calculated as the volume-average 
quantity. As the mass of dust in an object is proportional to pd 
(see Eqs. and Q), the assumption of porous grains can lead 
to considerably smaller mass estimates. Note that the opacities 
are larger for more absorbing carbon particles. A similar effect 
was noted by Quinten et al. (2002i who theoretically studied 
the wavelength dependence of extinction of different carbona- 
ceous particles. They also showed that the far IR extinction was 
larger for clusters of spheres and spheroids than for compact 
spheres. A very large enhancement of the submm opacities was 
found by Ossenkopf & Henning fl994J in the case of pure car- 
bon aggregates or carbon on silicate grains. 

Using Eq. (jfji and data from Table [5] shows how the parti- 
cle porosity and structure can influence estimates of dust mass 
in an object. The mass ratio estimates can be found in the 
Rayleigh-Jeans approximation 

Md(compact) KporousC-l) Td , porous 



Md(pOTOUS) Kcompact(/i) T d 

.compact 

With grain temperatures from Fig. [S] and the values of k for 
composite grains and EMT-Mie theory (the third column in 
Table we can find that the ratio Md(compact)/M d CP = 0.9) 
lies between ~ 1.3 and ~ 1.5. If the layered-sphere model 
is used (the forth column in Table [5} the ratio increases to 
3.8 - 4.3. This means that the calculated mass of an object can 
be reduced if compact grains are replaced by porous ones. 
The ratio of dust masses calculated for two grain models is 



Md(EMT-Mie) 



^lay sphere ,lay sphe: 



M d (layered sphere) k E mi-Uk{A) 7d,EMT-Mie 



2.8. 



The numerical value was obtained for particles consisting of 
AC1, astrosil and vacuum with f = 0.9 and the temperature 
ratio 7d,iaysphere/7d,EMT-Mie = 0.85 as discussed in Sect. 14. II If 
we consider particles of the same porosity but consisting of two 
materials, the ratio of masses will be even larger (3.1 for AC1 
and 4.3 for astrosil). Thus, one can overestimate the mass of an 



object by a factor of 3 or more if the EMT-Mie model is ap- 
plied, since real dust grains in molecular cloud cores should be 
very porous and should have non-Rayleigh inclusions. Another 
case when the effect can be important is in circumstellar discs, 
e.g. Takeushi et al. (2005 ) used k — 0.3 cm 2 /g at A = 1 mm for 
highly porous silicate grains, which is a good approximation 
only for particles with small size inclusions (see Table|5J. 



5. Concluding remarks 

We have considered how the porosity of composite cosmic 
dust grains can affect their optical properties important for in- 
terpretation of observations of interstellar, circumstellar and 
cometary dust. Two models of particle structure were used. 
Particles of the first kind had well-mixed inclusions small in 
comparison with the wavelength, while those of the second 
kind consisted of very thin, cyclically repeating layers. Earlier 
we showed that the optical properties of such layered particles 
are close to those of particles with small and large inclusions 
(see Voshchinnikov et al. I2005i . As effective medium theo- 
ries give reliable results for particles with small inclusions, two 
very different particle structure models can be simply realized 
and extensive computations can be performed. 

For both models, we studied how an increase of the volume 
fraction of vacuum could change the extinction efficiencies at 
different wavelengths, temperature of dust grains, profiles of 
the IR silicate bands and dust millimeter opacities. It is found 
that the models begin to differ essentially when the porosity 
exceeds ~ 0.5. This difference appears as lower temperatures 
(Sect. 14. it . shifted central peaks of the silicate bands ("Sect. l4.2t 
and larger millimeter opacities (Sect. I4.3i for layered particle 
model in comparison with that based on EMT calculations. The 
latter model also requires larger dust-phase abundances than 
the layered model (Sect. 13 .2t to produce the same interstellar 
extinction. 

The assumption that interstellar particles have only small 
size inclusions looks to some extent artificial (excluding, of 
course, the case of special laboratory samples). Therefore, we 
believe that the layered sphere model well describing light scat- 
tering by very porous quasispherical particles with inclusions 
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of different sizes should find wide applications in interpreta- 
tion of different phenomena. In particular, this model has good 
perspectives to explanation of the flat interstellar extinction ob- 
served in the near-IR part of spectrum (Sect. 13.31 and variations 
of the shape of the silicate feature detected in spectra of T Tau 
and Herbig Ae/Be stars (the results will be in a subsequent next 
paper). 
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